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ABSTRACT 

We study the spectral evolution on second and sub-second timescales in 1 1 long and 12 short 
Gamma Ray Bursts (GRBs) with peak flux >8.5 x 10~ 6 erg/cm 2 s (8 keV-35 MeV) detected 
by the Fermi satellite. The peak flux correlates with the time-averaged peak energy in both 
classes of bursts. The peak energy evolution, as a function of time, tracks the evolution of the 
flux on short timescales in both short and long GRBs. We do not find evidence of an hard- 
to-soft spectral evolution. While short GRBs have observed peak energies larger than few 
MeV during most of their evolution, long GRBs can start with a softer peak energy (of few 
hundreds keV) and become as hard as short ones (i.e. with Speak larger than few MeV) at the 
peak of their light curve. Six GRBs in our sample have a measured redshift. In these few cases 
we find that their correlations between the rest frame Speak and the luminosity L- lso are less 
scattered than their correlations in the observer frame between the peak energy E°^l k and the 
flux P. We find that the rest frame Speak of long bursts can be as high or even larger than that 
of short GRBs and that short and long GRBs follow the same S poa k(<) — Li SO (t) correlation, 
despite the fact that they likely have different progenitors. 

Key words: Gamma-ray: bursts — Radiation mechanisms: non thermal 



1 INTRODUCTION 

Since the discovery of the existence of the two classes of short and 
long GRBs, extensive studies of their prompt and afterglow prop- 
erties have been finalized to outline similarities and differences. 
Short GRBs have optical and X-ray afterglows similar to those of 
long GRBs (e.g. Gehrels et al. 2008, Nysewander, Fruchter & Pe'er 
2009) but, although still in a handful of cases, they are typically as- 
sociated with galaxies of almost all types (e.g. Berger et al. 2009) at 
preferentially lower redshifts with respect to those of long events. 
The leading idea is that short and long GRBs have different progen- 
itors although the mechanism producing their prompt and afterglow 
emission should be similar (e.g. see Lee & Ramirez Ruiz 2007 for 
a recent review). 

The comparison of the prompt 7-ray emission showed that 
short GRBs are harder than long ones, if the hardness ratio is used 
as a proxy of the spectral hardness (Kouveliotou et al. 1993). Con- 
sidering the time integrated spectra of samples of short and long 
GRBs in the BATSE population, Ghirlanda et al. (2009) found that 
short GRBs have a harder low energy spectral index while they 
have similar peak energy -E^lk compared to long bursts. Both in 
short and long GRBs, correlations are found between -Bp^ak an d 
the fluence (F) and between E°^l k and the peak flux (P) (Lloyd 
& Petrosian 1999, Ghirlanda et al. 2009, Nava et al. 2008, Nava 
et al. 2011b). In particular, while the E°^ k -F correlation is dif- 
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ferent for short and long GRBs (with short events having lower 
fluences for the same E°^l k of long bursts), they follow the very 
same -Epc ak — P correlation. Fluence and peak flux are the prop- 
erties commonly used to select samples of long and short bursts in 
comparative studies. Different selections, based on P or F, influ- 
ence the resulting distributions of -Epc a k of the selected samples 
and lead to different conclusions on the short-long -Epclk compari- 
son. Selection based on the peak flux P (Ghirlanda et al. 2009) re- 
sults in similar S°jJ k distributions for long and short GRBs, while 
a selection based on the fluence F leads to the conclusion that short 
bursts have on average larger peak energies (see Nava et al. 201 lb 
for a detailed discussion of this effect). 

Recently, Nava et al. (2011a) analysed the time integrated 
spectral properties of 438 GRBs detected by the Fermi satellite. 
Through this catalog, Nava et al. (2011b - Nllb hereafter) con- 
firmed previous results, based on the BATSE GRB population, of 
the existence of an -Epc ak — P ar, d an -^pcak — P correlation also in 
the Fermi sample and of the different spectral properties of short 
and long GRBs. Thanks to the wider energy range explored by 
the Gamma Burst Monitor (GBM, 8keV-35MeV) on board Fermi, 
Nllb also found that the peak energy -Epc a k of Fermi short bursts 
can reach very high values (confirming, with a larger population of 
short GRBs, what found for the three brightest short events detected 
by Fermi - Guiriec et al. 2010). Such high values of E^,% k are not 
so common in long bursts. However, these comparisons are typi- 
cally done in terms of observed peak energies. It is still possible 
that the different redshift distribution of the two classes of GRBs 
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(e.g. Guetta & Piran 2005) is at the origin of the difference in their 
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peak 



distributions. 



It was known that long BATSE GRBs show a strong spectral 
evolution (Ford et al. 1995; Ghirlanda, Celotti & Ghisellini 2002; 
Kaneko et al. 2006). This was confirmed also through time resolved 
spectral analysis of Swift bursts (Firmani et al. 2009). The study of 
the spectral evolution of Fermi long GRBs with measured redshifts 
(Ghirlanda et al. 2010, G10 hereafter) showed that a similar time 
resolved E°^l k (t) — P(t) correlation exists within these sources. 
Such a positive correlation means that the observed peak energy 
E°^l k tracks the flux P during the lifetime of the GRB. The knowl- 
edge of z for these events allowed to verify that a correlation in the 
rest frame similar to the Speak — Li so (so called "Yonetoku" corre- 
lation - originally discovered considering the GRB time-integrated 
spectra - Yonetoku et al. 2004) exists within individual bursts. A 
similar study on a sample of Fermi short GRBs (Ghirlanda et al. 
2011, Gil hereafter) confirmed that also in this class there exists 
a time resolved E°^l k (t) — P{t) correlation. These results suggest 
that a mechanism could be responsible for this internal (i.e. time 
dependent) correlation between the peak energy and the flux and 
that this mechanism is similar in long and short GRBs. 

However, the sample of G10 considered 12 long Fermi bursts 
with measured redshifts (until the end of July 2009) while the sam- 
ple of Gil considered short Fermi bursts selected on the basis of 
their peak flux and fluence (collected from the GCN circulars). 
Moreover, the analysis of G10 considered only the low energy Nal 
detectors (8 keV-1 MeV) of the GBM instrument, missing the pos- 
sibility to measure high E°^l k , while the spectral analysis of short 
GRBs of Gl 1 included also the BGO detectors that extend the spec- 
tral coverage to ~35 MeV. 

A direct comparison of the spectral evolution of short and long 
GRBs detected by Fermi requires a similar selection of the two 
samples (see e.g. Nava et al. 2011b for a discussion of the differ- 
ent selection criteria adopted in the literature to construct short and 
long GRB samples) and an homogeneous spectral analysis over the 
wider energy spectral range accessible by the GBM instrument, i.e. 
8 keV-35 MeV. 



2 THE SAMPLE 

The first compilation of the spectral properties of GRBs detected 
by Fermi was recently published in Nava et al. (201 la, Nl la here- 
after). They analyzed the public data of 438 GRBs detected by 
the Fermi GBM up to March 2010. For 318 cases they could reli- 
ably constrain the peak energy -E^k °f the V F V spectrum by ana- 
lyzing the time integrated spectrum between 8 keV and 35 MeV. 
Among these 318 GRBs, 274 are long and 44 are short. Long 
GRBs have a typical peak energy _E°^ k ~160 keV and low en- 
ergy photon spectral index a ~ —0.92. Short GRBs have harder 
peak energy (£pcak~490 keV) and harder low energy spectral in- 
dex (q ~ -0.50). 

Nlla also analyzed the spectrum corresponding to the peak 
flux of each GRB and in 235 (out of 438) cases they could constrain 
the peak energy E°^ k . Among these (shown in Fig.[TJi there are 31 
short (squares in Fig.QJ and 204 long (circles in Fig.[T} events. The 
peak spectrum of a burst typically corresponds to the time interval 
of the maximum of the flux in the light curve. If the spectral evo- 
lution of the GRB peak energy is correlated with the flux (as found 
independently for short and long Fermi bursts in G10 and Gl 1) then 
the selection on the peak flux ensures that we are considering the 
hardest short and long GRBs in the analyzed sample. Moreover, the 
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Figure 1. Observer frame peak energy versus peak flux for 235 Fermi 
GRBs analyzed by N10. Short GRBs (31 events) are shown with squares 
and long GRBs (204 events) with circles. The sample selected for this work 
comprises 11 long GRBs (filled red circles) and 12 short GRBs (filled blue 
square) with peak flux (8 keV-35 MeV) larger than 8.5xl0 _B erg/cm 2 s 
(shown by the vertical dashed line). 



typical selection of samples of long and short GRBs present in the 
literature adopts the peak count rate which is a rather poor proxy for 
the actual energy flux. The knowledge of the spectral shape for all 
bursts in the sample allows us to directly use the peak flux in physi- 
cal units (8 keV-35 MeV energy flux in erg/cm 2 /s) which accounts 
for the shape of the spectrum and for the peak energy and ensures 
to select at the same time the brightest and the hardest bursts both 
in the short and long population. 

We show in Fig. Q] the two selected samples of short (filled 
blue squares) and long (filled red circles) GRBs with P(8 keV-35 
MeV)^8.5xlO -6 erg/cm 2 s. The names of the selected bursts are 
reported in Tab. Q] together with their duration. The latter is taken 
from the table of Nlla. We note that this is not the Tgo parame- 
ter typically used to characterize the GRB duration, but the time 
interval over which the spectrum was accumulated and analyzed 
in Nlla. This time was selected by-eye from the light curve of 
each GRB as the time interval in which the signal was above the 
background. We caution that this criterion may result in a dura- 
tion which is smaller than the Tgo- However, all the short GRBs 
selected (except the short/hard burst 090510 detected by the Large 
Area Telescope on board Fermi - Abdo et al. 2009) have a duration 
much smaller than 2 seconds which is the empirical dividing line 
between short and long GRBs. 



3 SPECTRAL ANALYSIS 

The GBM (Meegan et al. 2009) comprises 12 thallium sodium io- 
dide [Nal(Tl)] and two bismuth germanate (BGO) scintillation de- 
tectors which cover the energy ranges ~8 keV-1 MeV and ~300 
keV-40 MeV, respectively. The GBM acquires different data types 
for the spectral analysis (Meegan et al. 2009) among which the 
"TTE" event data files containing individual photons with time and 
energy tags. Here we analyze TTE data for both short and long 
GRB because they allow to perform a time resolved spectral anal- 
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Table 1. Fermi GRBs selected for the time resolved spectral analysis (see 
Fig[TJ. Detectors analyzed: N=NaI and B=BGO. The last column shows the 
minimum integration time (i.e. maximum time resolution) adopted for the 
time resolved spectral analysis. 
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Figure 2. Peak energy versus flux for 749 time resolved spectra of 1 1 long 
GRBs (filled red dots) and 175 time resolved spectra of 12 short GRBs 
(filled blue squares). 



ysis on sub-second timescales for short GRBs. CSPEC data, with 
a maximum time resolution of 1.024s, are not suitable for time re- 
solved spectral analysis of short GRBs with typical duration <2 

s. 

For the time resolved spectral analysis we used the recently 
released software RMFn0(v33pr7). In order to model the back- 

2 http://fermi.gsfc.nasa.gov/ssc/data/analysis/user/ 
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Figure 3. Slope a of the correlation E°^ k (t) oc P(t) a for individual 
bursts. Filled red circles correspond to long GRBs, filled blue squares to 
short. The dot-dashed and dashed lines represent the average slopes for 
long and short GRBs, respectively. The solid filled and shaded regions are 
the lcr uncertainty on the average slopes. 



ground spectrum for the time resolved spectral analysis, we se- 
lected two time intervals before and after the burst. The sequence 
of background spectra in the two selected intervals were fitted with 
a first order polynomial to account for the possible time variation of 
the background spectrum. Then the background spectrum was ex- 
trapolated to the time intervals selected for the time resolved spec- 
troscopy of individual GRBs. 

We jointly fitted the spectra from at least 2 Nal and 1 BGO 
detector which had the largest illumination by the GRB. The in- 
clusion of the BGO data extends the spectral coverage of the Nal 
detectors from 1 MeV to ~35 MeV. Particularly bright short bursts 
(090227, 090228, 090328) were analyzed including more than 2 
Nal (and for 090510 both BGO detectors) as done by Guirec et al. 
(2010). We performed a time resolved spectral analysis by progres- 
sively reducing the time bin duration. For the brightest GRBs of 
our sample (i.e. GRB 090510 and GRB 090227 both of the short 
class) we performed a time resolved spectral analysis down to 4 
ms. In general the time resolution within a GRB is higher near the 
peak and lower along the tails of the burst. This was determined by 
the requirement that all the parameters of the fitted model to each 
time resolved spectrum could be constrained (i.e. their associated 
statistical errors determined with an accuracy < 100 %). In order 
to have a reasonable number of time resolved spectra and to make 
the plot of Fig. [2] readable, we limited the time resolution to 1024 
or 2048 ms for particularly long GRBs. In Tab.Q](last column) we 
report for each GRB analyzed the shortest time bin on which their 
spectrum was analyzed. 

The model adopted for the time resolved spectral analysis is a 
power law with an exponential cutoff whose free parameters are the 
low energy spectral index a, the peak energy -E^k (i- e the peak 
of the vF v spectrum), and the normalization. The choice of this 
model is motivated by the fact that at high energies the response of 
the BGO rapidly decreases for increasing energy so that it is hard, 
in single time resolved spectra, to constrain the possible presence 
of a power law component of e.g. the Band function (Band et al. 
1993) which is instead typically fitted to the time integrated spectra 
of long GRBs (see also Nl lb for a discussion on the model choice). 
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Figure 4. Left panel: observer frame peak energy versus peak flux for the analyzed GRBs. The 5 long and 1 short bursts with known redshifts are marked 
with different colors. Right panel: rest frame peak energy versus luminosity for bursts with measured redshifts. Solid and dashed lines are the -E pca k — Aso 
correlation and its 3<r scatter respectively (from Ghirlanda et al. 2010). 



The choice of the same model for the time resolved analysis of short 
and long GRB spectra ensures that if the model is introducing any 
systematic effect in the estimate of some of its spectral parameters, 
this is common to all the analyzed spectra. Note that if short and 
long GRBs spectra have intrinsically different shapes (e.g. a Band 
model for long and a cutoff-powerlaw for short) our choice would 
introduce a systematic effect, i.e. an overestimate of the peak en- 
ergy of long GRBs (e.g. Kaneko et al. 2006). Still it is hard, in time 
resolved spectra, to distinguish between a cutoff-powerlaw model 
and a Band function given the lower response of the instrument at 
high energies and the low fluxes of single time resolved spectra. 



4 RESULTS 

By considering the evolution of the peak energy -B^ak of the vF v 
spectrum and the energy flux P (integrated over the 8 keV-35 MeV 
energy range), we find that in most bursts there is a tracking pat- 
tern. This is shown by the strong correlation between E°^l k and 
P shown in Fig. [2] This correlation is shown for both short (filled 
blue squares) and long (filled red circles) GRBs and extends over 
four and three decades in P and -B^ak respectively. By consider- 
ing long and short GRBs separately, the Spearman correlation co- 
efficient is r=0.54 (with a chance probability of 10 -58 ) and r=0.74 
(with a chance probability of 10~ 31 ) for long and short GRBs, re- 
spectively. 

We also fitted the £°e ak (i) - P(t) correlation within the 
12 short and 11 long GRBs considered individually. The slopes 
a of the correlation E°^l k (t) oc P(t) a for individual bursts are 
shown in Fig. [3] with different symbols for long (filled red circles) 
and short (filled blue squares) GRBs. Fig. [3] also shows the slope 
(dashed and dot-dashed lines) obtained by averaging the slopes of 
short and long GRBs. The slopes of the E°^(t) - P(t) correla- 
tion in the sample of short GRBs are more scattered, with also one 
case (i.e. GRB 081209) having a very flat correlation. In general the 
correlation slopes of short and long GRBs are similar (as shown by 
their average values) and consistent with a typical slope of 0.5 also 
found within the sample of Fermi GRBs by considering the time in- 
tegrated spectra (Nllb). Moreover, this is also consistent with the 
slope of the "classical" Yonetoku correlation, i5 pca k — Li SO , (which 



is computed in the rest frame for bursts with measured redshifts) 
obtained considering the time integrated spectra of long and short 
GRBs (see e.g. G10 for a recent compilation of this correlation). 

Short and long GRBs follow the same correlation but they are 
slightly displaced in Fig. [2] with short GRBs having larger fluxes 
and peak energies with respect to long GRBs, similarly to what 
found with the time integrated spectral properties of Fermi and 
BATSE long and short GRBs (Nl lb). However, among the sample 
of long GRBs there is GRB 0809 16C whose peak energy E°^ k 
reaches values of few tens of MeV at the peak and becomes as hard 
as short GRBs (better shown in Fig. [4]- left panel). 

The quantities plotted in Fig.|2]are in the observer frame. Since 
the average redshift of long GRBs is larger than that of short ones, 
the displacement shown by Fig.|2]between short and long is prob- 
ably reduced in the rest frame. In our samples of bursts there are 6 
bursts with known z. Among these one is the short GRB 090510 at 
z=0.903. In Fig. [4] we highlight with colors and different symbols 
these bursts in the observer frame (left panel) and in the rest frame 
(right panel). We notice that when transforming from the observer 
frame E°^ v (t) - P(t) to the rest frame -E pcak (£) - L iso (t) the 
scatter of the correlations of these 6 bursts is evidently reduced. 
Moreover, in the rest frame (right panel in Fig|4]l the long GRB 
080916C (at z — 4.35) has the largest peak energies (i.e. larger 
than the short GRB 090510). 



5 DISCUSSION AND CONCLUSIONS 

Compared to previous analysis (G10, Gil and Guiriec et al. 2010) 
that considered separately short and long GRBs selected for their 
peak flux or fluence, here we ensure to have a complete sample 
down to our selection threshold. Moreover, we performed a homo- 
geneous time resolved analysis of these bursts over the 8 keV - 35 
MeV energy range by combining the Nal and BGO detectors of the 
GBM instrument on board Fermi. 

We find that both short and long GRBs follow a strong 
.Epcak (*) — P(t) correlation in the observer frame (Fig.[2j. Consid- 
ering single bursts, the slopes a of the E°^l k (t) oc P(t) a correla- 
tion are consistent with a a ~0.5 (with short GRBs having a larger 
scatter of the values of a than long ones, see Fig.[3j- 
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The existence of a E°^ k (t) — P(t) correlation in short and 
long GRBs indicates that the typical spectral evolution in both 
classes of events is tracking rather than hard-to-soft. This result 
is different from the dominating hard-to-soft spectral evolution 
found in a sample of BATSE bursts (Hakkila & Preece 2011) and 
in agreement with past studies based on the time resolved spectral 
modeling of bright BATSE GRBs (Ford et al. 1996, Ghiiianda et 
al. 2003). A possible reason could be that here we model the sin- 
gle time resolved spectra in each GRB with a spectral model (the 
cutoff-power law function) and use the resulting E°^l k as a proxy 
of the GRB hardness instead of the count hardness ratio as done in 
Hakkila & Preece (2011). 

Only 6 GRBs in our sample have measured reshifts and, 
among these, we have only one short GRB. Fig.|4]shows that in the 
Speak (t) — Liso(i) plane these bursts have a reduced scatter with 
respect to the E°^l k (t) — P{t) plane. Moreover, in the rest frame 
plane, long bursts can reach values of i? pca k as large as those of 
the only short one with known z, contrary to what occurs for the 
observed values (see e.g. Guirec et al. 2010). If these results, now 
based on a very limited number of events, would be confirmed by 
the addition of short and long GRBs with measured redshifts with 
a time resolved spectral analysis, it could indicate that short bursts 
appear harder (i.e. larger E°^l k ), than long GRBs, only because 
of their different average redshift. Interestingly, the 6 GRBs with 
measured redshift have individual time resolved -E poa k (t) — Li so (t) 
correlations which are all consistent with the -B pca k — Liso correla- 
tion (solid line in FigO found with the time integrated spectra of a 
much larger sample of GRBs. 

Several interpretations of the L pca k — Li so correlations have 
been proposed so far: (a) kinematic interpretations in which the 
link between _E poa k and Li so is established by the configuration of 
the emission region (Yamazaki et al. 2004; Nakamura 2000; Ku- 
mar & Piran 2000; Toma et al. 2005; Eichler & Levinson 2004); 
(b) radiative interpretations in which it is the emission mechanism 
of the prompt phase to link £v pea k and Liso, as in the case of a 
spectrum dominated by a thermal component (Ryde et al. 2006; 
Thompson, Meszaros & Rees 2007) or in the case of magnetic 
reconnection (Giannios & Spruit, 2007) or (c) both radiative and 
geometric effects (Lazzati et al. 2010) . All these interpretations 
however considered the time integrated L pca k — Li so correlation 
found in the population of long GRBs. Here we have shown that 
two new ingredients should be considered: (1) that a time resolved 
L" P e a k(£) — Liso(i) correlation exists within individual GRBs and 
that (2) this is similar in the population of long and short GRBs. 
Among the above mentioned models, Giannios & Spruit (2007) 
predict that a time resolved E peak (t) — Li so (t) correlation should 
exists within single GRBs and that more luminous GRBs should 
have a larger baryon loading. 

Our results indicate that both short and long GRB have in com- 
mon a similar physical mechanism that link the flux and the peak 
energy within individual bursts as a function of time. If short and 
long GRBs have a different progenitor (the merger of two compact 
objects and the death of a massive star respectively) the similar 
E peak (t) — Liso(i) correlation found in both populations should 
arise from the similar emission mechanisms (independent of the 
type of progenitor). We have shown that within individual GRBs 
(when considering time resolved spectra) there is strong correla- 
tion between the peak energy and the peak flux (or luminosity for 
the few events with measured redshifts) and this correlation has 
the same slope of the correlation defined among different GRBs 
(by considering their time integrated spectral properties - either the 
peak energy-peak flux correlation or the peak energy-luminosity 



Yonetoku correlation). The found E°^l k (t) — P{t) correlation can- 
not be due to the detector threshold, and its coincidence with the 
L p c a k — P correlation found among different GRBs suggests that 
the latter is not due to instrumental selection effects but should have 
instead a physical origin. 
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